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Summary 

Thrust-vectoring engine nozzles raise new pos- 
sibilities for controlling future generations of jet 
airplanes Aerodynamic control surfaces are conven- 
tionally 'used to generate the moments required to 
pitch, roll, and yaw airplanes. These moments be- 
come weak during low-speed flight because of low 
dynamic pressure. At high angles of attack, aerody- 
namic controls are corrupted by cross-axis coupling 
terms, which complicate the design of control sys- 
tems Thrust-vectoring controls, on the other hand, 
are especially effective when the dynamic pressure is 
relatively low. The moments generated by thrust- 
vectoring controls remain aligned with the axes of 
the airplane regardless of the angle of attack. The 
incorporation of thrust-vectoring controls into jet 
airplane designs promises to extend the low-speed 
region of the flight envelope and may influence in- 
flight maneuver capabilities and airfield operations 
In addition, thrust-vectoring controls provide redun- 
dancy for the aerodynamic controls, which is a sig- 
nificant advantage in the event of failures or dam- 
age. This report documents a study undertaken to 
explore the application of thrust-vectoring controls 
to the lateral and directional control of a jet fighter 
airplane. The result of the study is a preliminary 
feedback control system design for a specific airplane 
configuration. 

A mathematical model representative of a modem 
high-performance twm-jet fighter airplane including 
ailerons, rudder, and independent horizontal-tail sur- 
faces was formulated Idealized bidirectional thrust- 
vectoring nozzles capable of generating control mo- 
ments about the three airplane axes were added to 
this model. Linear perturbation equations for level, 
trimmed flight were calculated for angles of attack 
ranging from 0° to 90°. Scheduled gains between 
the control effectors and two pseudo control variables 
were found. These schedules were selected so that 
changes in a lateral pseudo control variable would 
cause coordinated aerodynamic and thrust-vectoring 
control motions which would influence the roll and 
spiral modes of the airplane with little disturbance 
of the Dutch roll mode In a similar fashion, a direc- 
tional pseudo control variable was chosen such that it 
would mainly influence the Dutch roll mode. A mod- 
ified system of perturbation equations was formed by 
combining these schedules with the original pertur- 
bation equations. The inputs of the modified system 
were the two pseudo control variables. Feedback and 
feedforward gains were obtained by solving a linear 
quadratic regulator problem using quadratic penal- 
ties on the differences between the airplane responses 
to lateral stick and rudder pedal movements and 


idealized responses. These gains were then scheduled 
as functions of angle of attack. 

Simulated responses to step pilot control inputs 
were stable and well behaved. For lateral stick 
deflections, peak stability axis roll rates were between 
1 25 and 1.60 rad/sec over a range of angle of attack 
of 10° to 70°. For rudder pedal deflections, the roll 
rates accompanying the sideslip responses could be 
arrested by small lateral stick motions. 

Introduction 

Lateral control and directional control of air- 
planes are traditionally treated as separate design 
problems. Rolling maneuvers (lateral control) are 
performed by lateral deflections of the control stick, 
which activate the ailerons to produce rolling mo- 
ments Sideslipping maneuvers (directional control) 
are performed by deflections of the rudder ped- 
als, which activate the rudder to produce yawing 
moments Acceptable control of most airplanes is 
achieved by mechanical systems which move the con- 
trol surfaces in proportion to the cockpit controls and 
by designing the airplane so that it is aerodynami- 
cally stable and well behaved. 

Modern high-performance airplanes which have 
wide flight envelopes are often difficult to handle with 
traditional control strategies in certain flight regimes 
Electronic systems are used in these airplanes to pro- 
duce artificial stability and improve handling qual- 
ities. Stability augmentation systems improve the 
damping of airplane motions, cause the airplane to 
respond better to the pilot’s control motions, make 
the expansion of flight envelopes possible, and per- 
mit innovations in airframe design. As these systems 
have evolved, they have remained segregated into lat- 
eral and directional systems even though there may 
be crossfeeds between them, for example, the aileron- 
rudder-interconnect used during flight at high angles 
of attack. 

Additional complications arise when new types of 
control effectors are added to yield unconventional 
airplane configurations. The ailerons may be aug- 
mented by spoilers and differential deflections of the 
horizontal tail, for example, and the rudder may 
be augmented by lateral thrust vectoring. Different 
control effectors are effective during different flight 
phases. For example, aerodynamic control surfaces 
can easily overpower thrust-vectored controls during 
high-speed flight but become inferior at low speeds. 

Controls which are effective for a given mode at 
one flight condition may affect a different mode at 
another flight condition. For example, the ailerons, 
which excite the roll mode during cruising flight, 
can cause directional oscillations during flight at 
high angles of attack because of excitation of the 



Dutch roll mode. Scheduling of the control devices 
is necessary to maintain desired airplane response 
to the pilot inputs throughout the flight envelope. 
The design of stability augmentation systems by 
the traditional axis-by-axis approach is becoming 
increasingly complicated because of (1) the widening 
range of flight conditions in which new airplanes must 
operate and (2) the large selection of control effectors 
available for new airframe designs. 

Modern multivariable controls methodology 
makes the simultaneous design of lateral and direc- 
tional controls possible. The stability and control 
characteristics of both axes are improved, and multi- 
ple redundant controls are utilized in the most ad- 
vantageous manner. Translation of the results of 
multivariable techniques into practical control sys- 
tem implementations is cumbersome The solutions 
are in the form of a matrix of gains, one gain for each 
state variable-control variable combination. Imple- 
mentation requires the measurement or estimation of 
all the state variables. Each gain may be expected to 
vary with changing flight condition. There is no dis- 
tinction between lateral and directional control sys- 
tems. Such a distinction is desirable because it per- 
mits analysis of system elements which can be related 
to recognizable items present in conventional control 
systems 

A preliminary design of a lateral-directional con- 
trol system for a fighter airplane utilizing thrust- 
vectoring engine nozzles is presented in this report. 
The incorporation of thrust-vectoring controls into 
jet airplane designs promises to extend the low-speed 
region of the flight envelope and may influence in- 
flight maneuver capabilities and airfield operations. 
The control system is specified for a wide angle-of- 
attack range, and the airplane has several redundant 
moment generators. This redundancy is a significant 
advantage in the event of failures or damage Such 
design problems are representative of the types that 
are currently challenging airplane designers. 

The present study maintains the distinction be- 
tween the lateral and directional control systems 
while utilizing multivariable design techniques to 
specify the system gains The control effectors are 
scheduled to form lateral and directional control 
channels by the application of the relative controls 
effectiveness technique of reference 1. Gains which 
specify the lateral and directional feedback signals 
are obtained by the solution of a linear quadratic reg- 
ulator (LQR) problem. Signals from onboard sensors 
are combined to form the feedback signals The feed- 
back gains are scheduled for operation over a wide 
range of flight conditions. Comparisons and rela- 
tive evaluations of the control system are beyond the 
scope of this report. 


Symbols 

A matrix of state variable coefficients 

a y measured lateral acceleration, m/sec 2 

ciy, eg lateral acceleration of center of gravity, 

m/sec 2 

B matrix of control variable coefficients 

b wing span, m 

C control-mixing matrix 

Cj) aerodynamic drag coefficient 

C e ff relative control effectiveness measure 

(see ref. 1) 

Ci aerodynamic lift coefficient 

Ci aerodynamic rolling-moment 

coefficient 

C m aerodynamic pitching-moment 

coefficient 

C n aerodynamic yawing-moment 

coefficient 

Cy aerodynamic side-force coefficient 

CCV control configured vehicle 

c mean aerodynamic chord, m 

D drag force, N 

F state variable feedback matrix 

f 6 measurement variable feedback matrix 

/ vector of coefficients relating state 

variables to lateral acceleration 

G matrix of pseudo control influences 

upon measurement vector 

g gravitational constant, 9.81 m/sec 2 

g vector of coefficients relating control 

variables to lateral acceleration 

H matrix of state variable influences 

upon measurement vector 

I identity matrix 

7xz product of inertia, kg-m 2 

I z moment of inertia, kg-m 2 

J performance index 

L lift force, N 
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LQR 

linear quadratic regulator 

Pest 

estimated rate of change of angle of 

lep 

distance of center of percussion from 


sideslip, rad/sec 


i center of gravity, m 

1 

flight path angle, deg 

k 

distance of tail from center of 

6a 

aileron deflection, deg (positive, right 


gravity, m 


aileron up) 

M 

pitching moment, N-m 

6d 

differential horizontal-tail deflection, 

m 

mass, kg 


deg (positive, right trailing edge up) 

P 

body axis roll rate, rad/sec 

6h 

horizontal-tail deflection, deg (positive, 

Pm 

stability axis roll rate model, rad/sec 


trailing edge down) 

Pa 

stability axis roll rate, rad/sec 

6r 

rudder deflection, deg (positive, 
trailing edge left) 

<? 

pitch rate, rad/sec 



■ dynamic pressure, Pa 

<5rp 

rudder pedal deflection, portion of 

Q 


full deflection (positive, left pedal 

r 

body axis yaw rate, rad/sec 


depressed) 

I'm 

stability axis yaw rate model, rad/sec 

6s 

control stick lateral deflection, portion 

s 

reference wing area, m 2 


of full deflection (positive for deflec- 
tion to right) 

s 

Laplace transform variable, sec -1 

6v, pitch 

pitching thrust vector, deg (positive, 

T 

thrust, N 

nozzles down) 

Ice 

turn coordination error, rad/sec 

^v,roll 

rolling thrust vector, deg (positive, 

t 

time, sec 

right nozzle up) 

u 

vector of control variables (denoted by 

6v, yaw 

yawing thrust vector, deg (positive, 


u on fig. 19) 

nozzles left) 

V 

true airspeed, m/sec 

Vl,V2 

Gaussian white noise variables 

Vi 

= max(50,P), m/sec 

6 

pitch angle, deg 

V 

vector of pseudo control variables 

A 

eigenvalue, sec -1 

w dir 

directional pseudo control variable 

t 

damping ratio 

«lat 

lateral pseudo control variable 

P 

air density, kg/m 3 

W 

weight, N 


bank angle, deg 

X 

vector of state variables (denoted by x 
■ on fig. 19) 

LJ n 

natural frequency, rad/sec 

Ym 

vector of measurements 

A dot 

(■) over a variable denotes differentiation 

yi 

difference between stability axis roll 

with respect to time. Derivatives of aerodynamic co- 

rate and roll rate model, rad/sec 

efficients with respect to angle of attack, sideslip an- 
gle, or any of the control deflections are denoted by 

V2 

difference between actual and desired 

subscripts as, for example, Ci a = dCi/da. Deriva- 


turn coordination error, rad/sec 

tives of aerodynamic coefficients with respect to pitch 

V3 

integrated error of lateral acceleration, 
m/sec 

rate are scaled by c/2V as, for example, C mq = 
(2V/c) dCm/dq. Derivatives of aerodynamic coeffi- 


cients with respect to roll rate and yaw rate are scaled 

a 

angle of attack, deg 

by b/2V 

as, for example, Ci = (2 V/b) dCi/dp. 

0 

angle of sideslip, deg 

The measure of filter gain (in decibels, dB) is cal- 
culated as 20 times the common logarithm of the 

Pest 

estimated angle of sideslip, rad 

gain. 
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Design Procedure Outline 

The steps used in the design procedure of devel- 
oping the lateral-directional control system for the 
fighter airplane utilizing thrust-vectoring engine noz- 
zles are briefly outlined in this section. More exhaus- 
tive elaboration is presented in following sections. 
The steps of the design procedure are 

1. Construct a mathematical dynamic model 

2. Find a set of trim conditions and perturbation 
equations. 

3. Form pseudo control variables. 

4. Solve a linear quadratic regulator problem 

5. Modify the feedback variable set. 

6 Delete ineffective gain elements. 

7. Schedule the gams. 

The unique feature of this process is the third 
step, in which five control effector variables are trans- 
formed into two pseudo control variables This step 
reduces the dimension of the gain matrix which must 
be processed in the remaining steps. 

The first step of the design procedure is to con- 
struct a mathematical model of the airplane dynam- 
ics. The aerodynamic coefficients used are simpli- 
fied functions of angle of attack. Generally, these 
were constructed of piecewise linear and low-order- 
polynomial segments selected to reasonably fit the 
available wind-tunnel data and to provide continu- 
ous derivatives where needed The aerodynamic co- 
efficients were defined for an angle-of-attack range of 
0° to 90°. 

Idealized bidirectional thrust-vectormg engine 
nozzles were appended to the airplane mathematical 
model to provide additional control moments. Pitch- 
ing moments were generated by simultaneously vec- 
toring the engine nozzles in the pitch plane of the 
airplane. Yawing moments were generated by simul- 
taneously vectoring the nozzles laterally. Rolling mo- 
ments were generated by differentially vectoring the 
nozzles in the pitch plane. This thrust-vectormg ca- 
pability permits controlled flight at high angles of 
attack and low airspeeds, at which the aerodynamic 
controls are ineffective 

The mathematical model was trimmed for hori- 
zontal flight at constant airspeed for angles of attack 
up to 90° (hovering flight). The horizontal tail was 
used to provide the necessary pitching moment until 
it reached its maximum deflection near the stall angle 
of attack At higher angles of attack, the thrust of the 
engines was vectored in the airplane pitch plane to 
supplement the horizontal tail. The dynamics of per- 
turbations about these trim conditions were found 
These linear perturbation dynamics were used in the 
development of the control system 


The next step was to determine schedules for 
the five lateral-directional controls (ailerons, rud- 
der, differential horizontal tail, and rolling-moment- 
producing and yawmg-moment-producmg thrust vec- 
tors). To accomplish this, the relative control effec- 
tiveness technique of reference 1 was applied to the 
perturbation equations. Two “pseudo” control vari- 
ables were found, with one concentrating its effec- 
tiveness on the roll and spiral modes and with the 
other primarily affecting the Dutch roll mode. Acti- 
vation of each of the pseudo control variables resulted 
in proportional deflections of the lateral-directional 
controls. The calculated ratios of the control deflec- 
tions to the pseudo control variables changed with 
flight condition These ratios were used as the basis 
for constructing schedules to interconnect the lateral- 
directional controls. Because of the uncertainty of 
having reliable air-data information available near or 
beyond stall, the schedules were frozen during flight 
in these regimes. The control schedules were com- 
bined with the linear perturbation models to form 
modified perturbation models whose inputs were the 
pseudo control variables The use of the pseudo con- 
trols minimized coupling of the modes of the airplane 
and reduced the number of control variables in the 
subsequent steps of the design procedure 

A linear quadratic regulator (LQR) problem was 
formulated to obtain feedback and feedforward gains 
which augmented the stability of the airplane and 
specified the signal paths from the cockpit controllers 
(stick and rudder pedals) and the state variables to 
the pseudo control variables. The modified pertur- 
bation equations were augmented by lateral and di- 
rectional pilot models. In each case the pilot was 
modelled by a Gaussian random variable driving a 
first-order lag. Two such first-order lags were in- 
cluded to provide idealized airplane lateral and di- 
rectional responses to the pilot model commands. 
The integral of the lateral acceleration was included 
to minimize long-term errors in turn coordination. 
The quadratic performance index for the LQR prob- 
lem included weights on the differences between the 
airplane responses and the idealized responses, the 
integral of the lateral acceleration, and the pseudo 
control variables. 

The gam matrix resulting from the solution of 
the LQR problem specified the commands to the 
pseudo control variables as linear combinations of the 
state variables of the augmented perturbation model. 
A transformed gain matrix was used to specify the 
pseudo control variables as linear combinations of 
the state variables, sensor variables, and other vari- 
ables which were functions of the state and sensor 
variables. The sideslip state variable was replaced 
by the lateral acceleration measured at the center of 
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percussion of the airplane relative to a pivot in the 
vicinity of the horizontal tail and engine nozzles. The 
bank angle was replaced by the calculated difference 
between the airplane turn rate and the ideal turn 
rate The ideal rate is a function of airplane attitude 
and velocity. 

The effectiveness of the elements of the modified 
gain matrix was examined at a number of angles 
of attack. Responses of the perturbation system to 
step deflections of the control stick and rudder pedals 
were obtained. Each element of the gain matrix was 
tested to determine if it could be removed without 
adversely affecting the airplane responses. Those 
elements which were determined not to be critical 
to the airplane’s stability and controllability by the 
pilot were eliminated 

The reduced gain matrices obtained in the pre- 
vious step were used to construct schedules The 
schedules were specified as piecewise linear functions 
of angle of attack. The airplane responses using the 
scheduled gains were checked by simulation at the an- 
gles of attack used in the previous step and at several 
intermediate angles. Because air data measurements 
may be unreliable near and beyond stall, the sched- 
ules were made constant in these flight regimes The 
signal formed by the difference between the airplane 
turn rate and the ideal turn rate based on airplane 
attitude and airspeed was found to be important at 
all the angles of attack examined However, because 
it depended on reliable measures of angle of attack, 
pitch angle, bank angle, and airspeed, it was consid- 
ered to be unavailable at high angles of attack. A 
substitute signal was formed by obtaining a deriva- 
tive of the lateral acceleration signal by means of a 
second-order differentiating filter. This signal was 
used at angles of attack above 20°. 

The design process outlined above could be ap- 
plied to the longitudinal perturbation system of the 
airplane. Given that a suitable control system is de- 
veloped for the longitudinal axis, the next step in 
the study would be to perform a piloted simulation 
study to evaluate the combined control system. Be- 
cause the addition of thrust-vectormg controls en- 
ables controlled flight into the poststall regime, new 
and unusual piloting techniques may surface to ex- 
ploit the new capability of the airplane. Modifica- 
tions of the system to deal with the nonlinear and 
coupled lateral-directional and longitudinal dynam- 
ics are to be expected. These topics are beyond the 
scope of this report and are being left for future 
research 

Design Details 

The previous section briefly outlined the salient 
features of a proposed procedure for the design of a 


lateral-directional stability augmentation and control 
system for a high-performance fighter airplane fitted 
with thrust-vectormg controls This section expands 
that discussion and presents each of the points men- 
tioned previously in greater detail. 

Mathematical Model 

A set of nonlinear differential equations was used 
to formulate a mathematical model of the six-degree- 
of-freedom dynamics of an airplane representative 
of a modern high-performance fighter. The aerody- 
namic data were approximated by manually fitting 
a series of linear and low-order-polynomial segments 
to wind-tunnel data. The data were available for an 
angle-of-attack range of 0° to 60°. The data were 
extrapolated to an angle of attack of 90°. The air- 
plane cannot sustain appreciable airspeed at such ex- 
treme angles of attack. The aerodynamic forces and 
moments would be very small in comparison with 
the forces and moments generated by the thrust- 
vectoring controls. For this reason and because the 
present study was not intended to deal with a spe- 
cific airplane, the data used were considered to be 
adequate for the present study. 

The lift and drag coefficients are depicted in 
figure 1. The lift curve is composed of a linear portion 
for angles of attack from 0° to 14°, a quadratic 
portion from 14° to 42°, and a linear portion from 
42° to 90°. This curve is smooth; therefore, the lift- 
curve slope is a continuous function of angle of attack. 
The drag coefficient is an exponential function of 
the lift coefficient at low angles of attack. Between 
24° and 60°, the drag coefficient is the product of 
the lift coefficient and the trigonometric tangent of 
the angle of attack. Above 60°, the drag coefficient 
is a quadratic function of angle of attack with a 
maximum value of 2.0 at 90°. 

The pitching-moment coefficient is depicted in 
figure 2. The pitching-moment curve is composed of 
a linear portion for angles of attack from 0° to 12°, 
a quadratic portion from 12° to 20°, and a linear 
portion from 20° to 90°. This curve is smooth; thus, 
the C ma derivative is a continuous function of angle 
of attack 

The effects of horizontal-tail deflections ( 6 h) upon 
the lift and pitching-moment coefficients are depicted 
in figure 3. The curves shown for Ci gh and C m $/, are 
composed of linear and quadratic functions of angle 
of attack. The pitch damping coefficient, Cmq, has a 
constant value of —10 6 (rad/sec) -1 

The lateral aerodynamic derivatives are depicted 
in figure 4, and the lateral control derivatives are 
depicted in figure 5. These derivatives are made up 
of piecewise linear functions of angle of attack with 
the exception of Ci ga , which includes a quadratic 
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segment. The sideslip derivatives C ^ and Cy^ 

are constant, with values of —0 00225 deg -1 and 
-0.0155 deg -1 , respectively. The side-force deriva- 
tives Cy and Cy r are 0. 

An idealized bidirectional thrust-vectoring nozzle 
was fitted to each of the airplane engines Each noz- 
zle was assumed to be capable of changing the di- 
rection of the thrust force by 15° in any direction. 
The operation of the two nozzles to produce moments 
about each axis of the airplane is depicted in figure 6 
Simultaneous deflections of the nozzles in the sym- 
metric plane of the airplane produce a pitching mo- 
ment (pitching thrust vector), differential deflections 
in the symmetric plane produce a rolling moment 
(rolling thrust vector), and simultaneous lateral de- 
flections produce a yawing moment (yawing thrust 
vector). The thrust force was assumed to pivot about 
the nozzle location with no change in magnitude. No 
consideration was made for variations in the aerody- 
namic coefficients caused by thrust vectoring. 

Perturbation Equations 

Trim solutions of the mathematical model were 
found for flight at constant speed and altitude over 
a range of angle of attack of 0° to 90°. The forces 
and moments acting on the airplane were balanced 
by deflections of the horizontal tail, pitching moment 
produced by thrust vectoring, and engine thrust, as 
depicted in figure 7. For low angles of attack (less 
than the stall angle of attack), the pitching moment 
required to trim the airplane was generated by the 
horizontal tail. At higher angles of attack for which 
the horizontal tail alone was unable to trim the 
airplane, it was supplemented by thrust vectoring. 
The thrust of the engines was calculated for flight at 
constant speed and altitude. Sea level standard air 
density was used for these calculations. 

The horizontal-tail deflection, 8 the total thrust 
of the engines, T, the thrust-vectoring angle, £„ iP itch> 
and the airspeed, V, of the trimmed airplane model 
are given in figure 8. At an angle of attack of 90°, 
the airplane is “standing on its tail” in a hovering 
condition. 

A system of linear perturbation equations for 
the lateral-directional dynamics of the airplane was 
calculated for angles of attack throughout the range 
of 0° to 90°. An example is given in table I. The state 
variables are the body axis roll and yaw rates, p and 
r, the bank angle, <f>, and the sideslip angle, /3. The 
control variables are aileron deflection, 8 a , rudder 
deflection, <5 r , differential horizontal-tail deflection, 
Sq, and the thrust-vectoring controls, 8 Vit - 0 u and 
6t>,yaw The lateral acceleration of the center of 
gravity (a y)C g) is included in the perturbation system. 


The open-loop eigenvalues of the perturbation 
system are plotted on the complex s-plane in fig- 
ure 9. For low angles of attack, the eigenvalues are 
configured m the traditional roll, spiral, Dutch roll 
arrangement. The roll and spiral modes are first or- 
der, with the eigenvalue of the spiral mode being near 
the origin and the eigenvalue of the roll mode being 
far from the origin on the negative real axis The 
complex eigenvalue shown represents the damped os- 
cillatory Dutch roll mode. At an angle of attack of 
approximately 20°, the roll and spiral modes are re- 
placed by a second-order mode, which is indicated 
on figure 9 by the real eigenvalues merging to form 
a complex pair This mode is commonly termed the 
“coupled roll-spiral” mode. At an angle of attack of 
40°, this mode becomes very lightly damped. The 
effects of this mode might be recognized by the pilot 
to be “wing rock.” 

It is difficult to describe the modes of the airplane 
in terms of the traditional modes (roll, spiral, etc.) 
when the angle of attack exceeds 20°. In this report, 
the mode corresponding to the complex eigenvalues 
which lie on the locus which originates with the 
merger of the roll and spiral eigenvalues is termed 
the “roll-spiral” mode The other second-order mode 
is termed the “Dutch roll” mode. 

The magnitudes of the lateral-directional eigen- 
values and the damping ratios of the second-order 
modes are plotted as functions of angle of attack in 
figure 10. For complex modes, the magnitudes shown 
are the undamped natural frequencies of the modes 
For real modes, the magnitude is the inverse of the 
time constant of the modes. 

The natural frequency of the Dutch roll mode is 
less than 1.5 rad/sec for angles of attack greater than 
15° and less than 1.0 rad/sec for angles of attack 
near 25° and greater than 77°. The damping ratio 
varies between 0.25 and 0 85 with the greater values 
occurring between 25° and 45° The maximum 
damping factor (£w n ) of 1.22 sec -1 occurs at an angle 
of attack of 40°. The damping factor is less than 
0.8 sec -1 for angles of attack between 5° and 26° 
and greater than 45°. 

For angles of attack less than 18°, the roll and 
spiral modes are distinct first-order modes. The 
time constant (|A| -1 ) of the roll mode increases with 
increasing angle of attack, whereas the time constant 
of the sprial mode decreases. Both time constants 
approach 2 sec at an angle of attack of approximately 
18°. For angles of attack greater than 18°, the roll 
and spiral modes are replaced by a second-order roll- 
spiral mode, which has a minimum damping ratio of 
0 094 at 40°. 
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Pseudo Control Variables 

The lateral-directional perturbation equations ob- 
tained in the preceding section include the effects of 
five independent control variables. The pilot can only 
be expected to be able to manipulate a maximum of 
two manual controllers for lateral-directional control 
of the airplane. These are typically lateral control 
stick and rudder pedal movements. Several sensed 
variables must be fed into the flight control comput- 
ers to stabilize the airplane and to modify its re- 
sponses in order to meet mission requirements and 
have acceptable handling qualities The combination 
of a large number of inputs and a large number of 
controlled outputs results in a sizable quantity of sig- 
nal paths which could be implemented through the 
control computers. In addition, the variation of the 
effectiveness of the aerodynamic and thrust-vectoring 
controls with changing flight condition makes the in- 
corporation of gain-scheduling algorithms into the 
flight control computers necessary. The resultant 
complexity of the final system could easily become 
quite unwieldy. 

The approach taken m the present study is to sys- 
tematically interconnect the controls, so that their 
activity is based upon two “pseudo” control vari- 
ables within the flight control computers The re- 
lationships between the controls and the pseudo con- 
trols are specified by schedules which are dependent 
upon flight condition. The schedules are developed 
from calculations based upon the effectiveness that 
each control has on the lateral-directional modes of 
the airplane. Incorporation of the schedules into the 
perturbation equations results in a modified system 
of perturbation equations. The modified system has 
two pseudo control variables for inputs, whereas the 
original system had five control inputs. The control 
interconnect schedules result m coordination among 
the original control variables. 

Before specifying the control interconnect sched- 
ules, it would be informative to examine the relative 
effectiveness of the five controls The relative control 
effectiveness measure of reference 1 was applied to 
the perturbation dynamic equations. This measure 
provides some insight into how the effectiveness of 
the different controls varies with changing flight con- 
dition. A control which has a large effect on a given 
mode is a logical candidate for use in (1) exciting that 
mode in response to pilot inputs and (2) modifying 
the frequency and damping of the mode by means of 
feedback loops. However, if the control must be used 
on a second mode, such usage will likely disturb the 
first mode. 

The relative effectiveness of each control upon 
the Dutch roll mode is presented in figure 11 as 


a function of angle of attack. The significance of 
this figure is that it quantitatively represents the 
ability of each control, relative to the total set of 
controls, to affect the Dutch roll mode. For angles of 
attack less than 16.5°, the Dutch roll mode is more 
strongly affected by deflections of the rudder than 
it is by any one of the other controls. For angles 
of attack greater than 16.5°, yawing thrust vector 
is the most effective. The rudder loses effectiveness 
at high angles of attack because of a small C ngr 
derivative and low dynamic pressure. The thrust- 
vectoring controls become powerful at high angles of 
attack because of the high engine thrust level At 
an angle of attack of 16.5°, the ailerons, rudder, and 
yawing thrust vector are equally effective upon the 
Dutch roll mode because the ailerons have an adverse 
yaw characteristic and because the Dutch roll mode 
involves large rolling motions about the longitudinal 
body axis of the airplane. The relative effectiveness 
of the rolling thrust vector reaches a level of 0 6 at 
an angle of attack of 70° (not shown) but remains 
less effective than the yawing thrust vector. At such 
an extreme angle of attack, the Dutch roll mode is 
characterized mainly by production of sideslip angle 
when the airplane rotates about its longitudinal axis. 
The rolling thrust vector is weak because of the close 
spacing of the engines, whereas the yawing thrust 
vector is much stronger because of the long moment 
arm between the airplane center of gravity and the 
engine nozzles. The strong yawing moments couple 
through the product of inertia, IxZi to excite the 
Dutch roll mode to a greater extent than the rolling 
moments, which act more directly on the mode. 

The relative effectiveness of the controls upon the 
roll and spiral modes is depicted in figure 12. At low 
angles of attack, aileron control is the most effective 
because it can produce the largest rolling moment. 
At higher angles of attack (up to approximately 16°), 
the rudder is the most effective because yawing mo- 
ments are required to rotate the airplane about the 
longitudinal stability axis (the wind vector) and be- 
cause the ailerons possess an adverse yaw characteris- 
tic At angles of attack greater than 16°, the yawing 
moment produced by the thrust-vectoring control is 
most effective in exciting the roll and spiral modes. 

The above discussion implies that the flight con- 
trol system should use the aerodynamic control sur- 
faces during high-speed flight at low angles of attack 
and the thrust-vectoring controls during low-speed 
flight at high angles of attack. The controls need 
to be blended, so that the most effective combina- 
tion is always in use. In addition, it is desirable that 
each control channel (lateral and directional) create 
a minimum disturbance on the other channel. To 
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accomplish these objectives, the control-mixing tech- 
nique of reference 1 has been used 

The control-mixing technique is a method for 
calculating proportional deflections of the controls, 
so that selected modes are affected to the max- 
imum extent possible while maintaining a mini- 
mum disturbance of the modes which are to re- 
main decoupled. This trade-off is accomplished by 
means of an algebraic optimization problem whose 
solution is the desired control combination. De- 
tails about the technique may be found in refer- 
ence 1. 

The control-mixing technique was applied to the 
perturbation system to find two sets of control in- 
terconnects The first set has maximum effect upon 
lateral airplane motions with minimum effect upon 
directional motions The second set has maximum 
effect upon directional airplane motions with mini- 
mum effect upon lateral motions The control inter- 
connects were calculated for the airplane with and 
without thrust-vectoring controls. 

The control interconnects calculated for lateral 
control (maximum excitation of the roll and spiral 
modes with minimum excitation of the Dutch roll 
mode) using only aerodynamic controls are presented 
m figure 13. At low angles of attack, lateral control 
is accomplished by simultaneous deflections of the 
ailerons and differential tail to produce rolling mo- 
ments At higher angles of attack, the rudder is used 
to produce the yawing moments required to roll the 
airplane about the longitudinal stability axis The 
ailerons reverse sign at angles of attack above 20°, 
so that their adverse yaw characteristic is utilized in 
aiding the rudder in producing yawing moments. In 
this region, the differential tail is producing rolling 
moments to cancel the rolling moments produced by 
the ailerons 

The control interconnects calculated for lateral 
control using thrust-vectoring controls in addition 
to aerodynamic controls are presented in figure 14. 
Below an angle of attack of 15°, the interconnects 
are similar to those presented in figure 13. The 
aerodynamic controls are much more effective than 
the thrust-vectoring controls when the airplane is 
trimmed in this region because of the high dynamic 
pressure The yawing thrust vector, £u,yaw, becomes 
effective between angles of attack of 10° and 20°. 
Above 20°, it overpowers the aerodynamic controls 
For this reason, the deflections of the aerodynamic 
controls are very small, whereas 5 v ,yaw operates at 
full authority. 

The calculated interconnects, if implemented in 
the airplane flight control computers, result in 
optimal decoupling of the airplane modes The cal- 
culated interconnects are useful as guides in speci- 
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fying practical control interconnect schedules. Inter- 
connect schedules were developed for lateral airplane 
control. These schedules, presented in figure 15, 
were manually fitted to the curves of figures 13 and 
14. The schedules for the aerodynamic controls were 
made to be functions of angle of attack and are based 
upon the curves of figure 13. The curves of figure 14 
were not used here, so that the aerodynamic controls 
would remain active at high angles of attack. A loss 
of thrust will not necessarily eliminate all possibil- 
ity of controlling the airplane. The thrust-vectoring 
controls were scheduled as functions of effective air- 
speed to approximate the curves of figure 14. Air- 
speed was used as the scheduling quantity instead 
of angle of attack, so that the airplane would re- 
tain effective controls if it should acquire a low angle 
of attack while at low airspeed during maneuvering 
flight. Alternatively, the thrust-vectoring controls 
could have been scheduled as functions of dynamic 
pressure. 

The interconnects calculated for directional con- 
trol (maximum excitation of the Dutch roll mode 
with minimum disturbance of the roll and spiral 
modes) using the aerodynamic controls are presented 
in figure 16 At low angles of attack, directional 
control is accomplished by the rudder. Between 10° 
and 20°, a combination of the rudder, ailerons, and 
differential tail is used. The deflections are in the 
crossed-control sense: rudder for nose-left moment, 
and ailerons and differential tail for right-wing-down 
moment. For angles of attack greater than 20°, si- 
multaneous aileron and differential tail deflections 
are used to produce a rolling moment about the air- 
plane longitudinal axis 

The interconnects calculated for directional con- 
trol using the thrust- vectoring controls in addition to 
the aerodynamic controls are presented in figure 17. 
These interconnects are very similar to those pre- 
sented in figure 16 The utilization of the rolling 
thrust vector gradually increases at angles of attack 
from 10° to 60°. This control is limited in its abil- 
ity to generate rolling moments because of the close 
spacing of the engines. The horizontal tail retains the 
ability to generate rolling moments at high angles of 
attack. Therefore, the transition between aerody- 
namic and thrust-vectoring controls is spread over a 
wide range of angle of attack. 

Control interconnect schedules were developed for 
directional airplane control. These schedules, pre- 
sented in figure 18, were manually fitted to the curves 
of figures 16 and 17. The aerodynamic controls 
were scheduled as functions of angle of attack and 
are based on the curves of figure 16 The thrust- 
vectormg controls were scheduled as functions of ef- 
fective airspeed and are based on the curves of fig- 



ure 17 The schedule for the rolling-moment thrust 
vector, (5„ r oii, was adjusted so that full authority was 
used for airspeeds ranging from zero to just greater 
than the stall speed. 

The schedules for the aerodynamic controls given 
in figures 15 and 18 require angle-of-attack sensor 
data ranging from zero to a value less than the stall 
angle of attack. At higher angles of attack, a stall 
sensor could be used. The schedules for the thrust- 
vectormg controls require airspeed sensor data for 
airspeeds above the stall speed. For slower speeds, 
accurate data are not required. 

This section has discussed the development of 
control interconnect schedules which specify how the 
five airplane controls need to be coordinated to result 
in decoupled (as much as possible) lateral and direc- 
tional control Two “pseudo” control variables can 
now be defined within the flight control computers 
The first pseudo control variable, acts through 
the schedules of figure 15 to deflect the controls in 
the proportions required to primarily excite the roll 
and spiral modes In the examples of reference 1, this 
pseudo control variable was treated as if it were a lat- 
eral deflection of the control stick In these examples, 
the airplane model responded to step inputs with a 
smooth buildup of bank angle and exhibited no ev- 
idence of adverse yaw motions The second pseudo 
control variable, v,j, r , acts through the schedules of 
figure 18 to primarily excite the Dutch roll mode. In 
reference 1, this pseudo control variable was treated 
as if it were a deflection of the rudder pedals. The 
airplane model responded to step rudder pedal inputs 
with improved sideslip behavior accompanied by less 
rolling motion compared with the response to rudder 
surface deflections alone. 

The perturbation system equations were com- 
bined with the control interconnect schedules to pro- 
duce the modified perturbation system shown in fig- 
ure 19. A numerical example is given in table II. The 
pseudo control variables are the control variables of 
the modified system. The control system design de- 
scribed in the remainder of this report is based upon 
the modified perturbation system. 

Linear Quadratic Regulator 

The airplane mathematical model used in the 
present study has poor stability at high angles of 
attack The airplane motions must be measured, 
and the controls must be driven by these measure- 
ments to improve stability. The design of a stability 
augmentation system involves the determination of 
the sensor complement which can make the required 
measurements and the development of control laws 
which transform the sensor signals into control vari- 


able commands Also, the control laws must trans- 
late the motions of the pilot’s controllers into control 
activity which results in acceptable airplane behav- 
ior. Although the results of reference 1 indicate that 
the control stick should be connected to V] at and the 
rudder pedals should be connected to w<j ir , this ar- 
rangement may not be appropriate in general In 
order to determine suitable gams to direct the sen- 
sor and pilot controller signals to the pseudo control 
variables, the linear quadratic regulator (LQR) tech- 
niques of reference 2 were applied to the modified 
perturbation system. This section describes the for- 
mulation of the LQR problem used in the present 
study. Following sections describe the process of 
making a more practical implementation of the con- 
trol laws. 

The LQR technique finds a matrix of feedback 
gains which specify the control variables of a linear 
system as linear combinations of its state variables 
The gains are calculated to minimize a quadratic 
function of system variables, which may include state 
variables, output variables, control variables, and 
combinations of these. The quadratic performance 
index is constructed to effect a trade-off between the 
magnitudes of undesirable system responses and the 
control variables. The result of the LQR technique 
is a gain matrix which makes the closed-loop system 
stable and which regulates the system to continually 
seek the nominal operating condition. 

The objectives of the LQR design for the present 
study are (1) to improve the lateral-directional sta- 
bility of the airplane, (2) to provide acceptable con- 
trol characteristics by using lateral control stick and 
rudder pedal motions, and (3) to coordinate turn- 
ing maneuvers. The third objective may be relaxed 
during operation at extremely high angles of attack, 
at which the airplane is not likely to be operated 
for sustained periods of time. The formulation of 
the LQR problem used to meet these objectives is 
described in this section. This approach closely par- 
allels the formulation used in the development of the 
control configured vehicle (CCV) system reported in 
reference 3 

A number of auxiliary equations were appended 
to the modified perturbation system. The added 
equations were lateral and directional pilot command 
models, lateral and directional reference response 
models, differences between the airplane motions and 
the reference models, and an integration over time of 
the lateral acceleration. These auxiliary equations 
are schematically depicted in figure 20. 

The pilot command models are first-order systems 
which are driven by Gaussian white noise. The 
break frequency of each is 1 rad/sec These models 
represent lateral control stick, 6 S , and rudder pedal, 
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<5 rp , motions which are limited to a bandwidth of 
1 rad/sec. Since the pilot controller motions are state 
variables, the LQR design technique calculates the 
pilot controller gains. 

The reference response models are first-order sys- 
tems which are driven by the pilot command models. 
The break frequency of each is 5 rad/sec. These mod- 
els respond to pilot commands to generate reference 
signals against which the airplane responses can be 
compared 

The first reference response model, p m , is the 
ideal airplane roll rate about the stability axis. The 
actual stability axis roll rate, p a , is a function of the 
body axis roll rate, p, the body axis yaw rate, r, and 
the angle of attack, a. The difference between the 
actual roll rate and the ideal roll rate is calculated 
by 

2/1 =Ps ~ Pm = pcos a + r sin a- p m (1) 

During coordinated turns, the turn rate is a func- 
tion of the component of the lift vector which is 
directed toward the center of the turn. This force 
is balanced by the centrifugal force of the airplane 
following a curved trajectory. This balance is not 
achieved during uncoordinated turns because of side 
forces, principally caused by sideslip, acting on the 
airplane The difference between the actual turn rate 
and the turn rate expected for a coordinated turn can 
be defined as the “turn coordination error,” T ce , as 
follows: 

T ce = y4>cos9 — (rcosa — psina) (2) 

where 

0 pitch angle 

(p bank angle 

V velocity 

g gravitational acceleration constant 

A variable which is interpreted as the rate of change 
of sideslip angle is defined in reference 3. This 
variable and T ce are calculated by identical formulas. 

For the present study, the response of the direc- 
tional reference response model, r m , to rudder pedal 
inputs was considered to be the desired response of 
the turn coordination error. The difference between 
the actual and the desired turn coordination error is 
calculated by 

2/2 = T ce - r m = ^tpcosd- (rcosa -psina) - r m 

( 3 ) 

The integral over time of the lateral acceleration 
of the airplane is a measure of the long-term error of 


turn coordination. If a directional control is applied 
and the wings are held level, the force produced by 
sideslip causes the airplane to execute a flat turn. 
The lateral acceleration and the turn rate are related 
through velocity. The difference between the lateral 
acceleration, a y , and the “ideal” acceleration, Vr m , 
is integrated over time to yield the desired measure 
as follows. 

2/3 = J{ay-Vrm)dt (4) 

The value of the lateral acceleration used here is 
taken at the location of the lateral accelerometer, 
which is discussed in a later section of this report. 

The performance index for the LQR problem is 
the integral over time of a weighted sum of squares 
of (1) the differences between the airplane responses 
and the idealized responses (yi, y 2 , and y%) and 
(2) the pseudo control variables, as follows: 

J = J (50 y\ + 2500 y\ + O.lyf + SOOwg* + MOt&Jdt 

( 5 ) 

The weights used in equation (5) are roughly equiv- 
alent to the weights used in reference 3. The LQR 
techniques of reference 2 were used to find the “feed- 
back” matrix of gams which minimize the perfor- 
mance index subj'ect to the dynamics of the modi- 
fied perturbation system and the auxiliary equations, 
which can be represented as 

x = Ax + Bv (6) 

where 

A matrix of state variable coefficients 

B matrix of control variable coefficients 

v pseudo control vector consisting of iq at 
and u dir 

x state vector consisting of (1) airplane 
state variables (p, 0, r, /?), (2) pilot 
command model states (<5 3 , <5rp), 

(3) reference model states ( Pm,r m ), 
and (4) integrated error of lateral 
acceleration (ys) 

An example solution of the LQR problem is given 
in table III. The solution is a matrix of gains which 
specify the pseudo control variables (uj at and v dir ) as 
linear combinations of all the state variables. Since 
the state variables include the pilot command mod- 
els and reference state models, the matrix includes 
“feedforward” gains for the pilot controllers. 
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Feedback Variable Change 

The sensor package carried onboard the airplane 
is assumed to include rate gyros to measure the body 
axis roll rate (p) and the yaw rate (r), a linear ac- 
celerometer oriented to measure the lateral accelera- 
tion {ay), and an attitude gyro to measure the bank 
angle (0). The lateral accelerometer is positioned 
forward of the center of gravity at the “center of per- 
cussion” of the airplane with respect to forces applied 
at the tail Considering only the mass and inertia of 
the airplane, application of control forces causes it to 
rotate about the center of percussion. The location 
of the center of percussion is the distance l cp forward 
of the center of gravity The formula for calculating 
lcp is 


where 

m mass of airplane 

Iz moment of inertia about vertical axis 

It distance between center of gravity and 

tail controls 

The accelerometer, with this positioning, princi- 
pally reacts to aerodynamic side forces (in this ex- 
ample, caused by sideslip acting through the aero- 
dynamic derivative Cy - 0 ) and has a minimum direct 
reaction to control motions. This location was chosen 
to avoid possible stability problems caused by a loop 
from control deflection to side force to accelerometer 
response. 

A measurement vector (y m ) was constructed of 
the known state variables ( 6 S , S ip , p m , r m , and 
1/3), the sensed angular rates (p and r), the sensed 
lateral acceleration {a y ), and the calculated value of 
the turn coordination error (T ce ). This vector is a 
linear function of the state vector of the perturbation 
system (x) and the pseudo control vector (v) and has 
the following form: 

y m = Hx + Gv (8) 

The solution of the LQR problem found in the pre- 
ceding section has the following form: 

v = Fx (9) 

Combining equations (8) and (9) yields 

v = F 6 y m (10) 

where 

F 6 = (1 + FH -1 G) -1 FH (11) 


The matrix (F&) is the feedback matrix, which spec- 
ifies the pseudo control variables as linear combina- 
tions of the elements of the measurement vector. Ex- 
amples are presented in table IV. The values of the 
elements of this matrix are plotted in figure 21. The 
data shown are the gains on the measurement vector 
which cause the controls to act in a manner which 
satisfies the LQR problem. 

The next step of the design process is to identify 
those elements of the gain matrix which are necessary 
for the satisfactory operation of the airplane controls. 
By the elimination of ineffectual gain elements, the 
mechanization of the control laws in the flight control 
computers can be simplified. 

Gain Element Reduction 

The individual gam elements of the feedback ma- 
trix were deleted or retained based upon the effects of 
their presence on simulated transient responses. The 
modified perturbation system was used to obtain the 
airplane responses to step lateral stick and rudder 
pedal movements. Reduced gain matrices were ob- 
tained at angles of attack of 10°, 20°, 30°, and 60°. 
This section describes the process of gain element re- 
duction and discusses the reduced-gain matrices. 

Simulated responses to step lateral stick and rud- 
der pedal movements using the full feedback matrices 
were obtained. When the stick was held off center, 
the airplane tended to return to a wings-level atti- 
tude after a long period of time This was caused by 
the feedback of 7/3 (the integral of lateral accelera- 
tion). Since the perturbation model did not include 
increasing pitching moment, which is required for 
turn coordination, the airplane model lost altitude 
when banked As a result, the simulation showed a 
buildup of sideslip and, consequently, lateral acceler- 
ation. The control system acted to drive the lateral 
acceleration to zero by causing the airplane to return 
to level flight Since this behavior is unacceptable in 
maneuvering flight, the feedback gains for 2/3 were 
deleted. These gains could be reintroduced into the 
control system as part of a more comprehensive sys- 
tem design. For the present study, these gains were 
not considered further. 

Simulated responses to step lateral stick and rud- 
der pedal movements were obtained with the con- 
troller crossfeeds (<5 S to t;j, r and <5 rp to i>j at ) and the 
reference model crossfeeds ( p m to w ( j ir and r m to U] at ) 
deleted. The elimination of these crossfeeds caused 
no significant deterioration of the airplane responses. 
These crossfeeds were not considered further in the 
present study. 

The elements of the feedback gam matrix which 
operate on signals derived from airplane motion sen- 
sors (p, T ce , r, and a y ) were examined next. Table IV 
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presents the measurement feedback gain matrices for 
the four angles of attack presently being considered 
The gain elements which are retained after the elimi- 
nation process are given in boldface type. In general, 
as the angle of attack is increased, more feedback 
loops must be utilized to control the airplane. 

For an angle of attack of 10°, most of the feedback 
gains were eliminated because either they had little 
discernible effect on the airplane responses or the 
effect was adverse. The only feedback gam element 
that was retained closed the loop from T ce to v<j ir 
This gain improves Dutch roll mode damping and 
causes a coordinating rudder deflection, which is 
desirable in turning flight. 

For an angle of attack of 20°, additional gams 
were used to close the loops from the yaw rate (r) and 
sensed lateral acceleration ( a y ) to the lateral pseudo 
control variable (uj at ) Both feedbacks improved the 
Dutch roll damping, with the a y loop having the 
greater effect. For a step rudder pedal deflection to 
the left, the airplane transitions smoothly to a steady 
right sideslip (nose-left) with a small right-wing- 
down bank angle. For a step lateral stick deflection, 
the airplane rolls smoothly to a constant bank angle 
with a time constant of approximately 4 sec. The 
airplane does not continue to roll because of the 
presence of the coupled roll-spiral mode and the yaw 
rate feedback The pilot must hold a stick deflection 
toward the turn in order to maintain a constant bank 
angle. The airplane returns to a wings-level attitude 
when the stick is released. Deletion of the yaw rate 
feedback results in a sustained rolling motion because 
the lateral acceleration feedback drives the controls 
to reinforce the pilot’s roll command This may be 
objectionable because the pilot must hold a stick 
deflection away from the turn in order to maintain 
a constant bank angle. Also, the control power 
available to right the airplane would be limited. A 
compromise value of 0 3 for the yaw rate to gain 
was selected This is approximately one half of the 
value given by the LQR problem. 

For an angle of attack of 30°, an additional gain 
element was used to close the feedback loop from 
the yaw rate (r) to the directional pseudo control 
variable (v<nr)- For a step rudder pedal deflection 
to the left, the airplane transitions smoothly to a 
steady right sideslip with a small right-wing-down 
bank angle For a step lateral stick movement, the 
airplane smoothly rolls to a constant bank angle with 
a time constant of approximately 25 seconds. The 
inclusion of the yaw rate to u^ir loop causes a rolling 
moment which lifts the wing on the side toward the 
turn. As a result, the sideslip angle and the lateral 
pseudo control variable (vi at ) are reduced. 


The implementation of the turn coordination er- 
ror (T ce ), used in formulating the measurement feed- 
back matrix in the preceding section, requires mea- 
surements of angle of attack, pitch angle, and bank 
angle. Obtaining these measurements reliably at 
high angles of attack, particularly near or beyond 
stall, may be impractical Also, operation near a 
pitch angle of 90° is likely to cause trouble because 
of the singularity existing in the Euler angles at this 
attitude. In order to avoid these problems, a substi- 
tute signal for T ce was constructed. 

As mentioned in the preceding section, one in- 
terpretation of the meaning of T ce is that it is an 
estimate of the time rate of change of the sideslip 
angle. The sideslip angle can be estimated from the 
sensed lateral acceleration of the airplane. For angles 
of attack of 30° and greater, the turn coordination 
error (T ce ) was replaced by an estimated value of the 
rate of change of the sideslip angle (/? e st) by using 
the following formula 

a _ s m 

(! + !)(! + ! )*sc Y0 ay (12) 

The first term of the formula is a second-order fil- 
ter which operates as a differentiator on signal com- 
ponents whose frequencies are less than 3 rad/sec. 
Higher frequency components are attenuated by the 
action of the second-order denominator, reducing the 
sensitivity of the differentiating filter to noise on 
the accelerometer signal but introducing phase lag, 
which may cause stability problems. Therefore, the 
use of this filter is limited to flight at low airspeeds, at 
which the closed-loop frequencies are low. The break 
frequencies of the filter denominator were selected 
from observations of their effects on the simulated 
airplane responses. 

The second term of equation (12) is the ratio 
between lateral acceleration and the sideslip angle 
based upon the value of the side-force derivative 
(CyA With the location of the lateral accelerometer 
at the center of percussion (discussed in the preceding 
section), the forces produced by the rudder and 
yawing thrust vector are not sensed by it. This 
location prevents objectionable coupling of the Dutch 
roll and roll-spiral modes. 

For an angle of attack of 60°, the solution of the 
LQR problem resulted in no feedback of yaw rate 
(r) to the directional pseudo control variable (udir)- 
However, a value of 0 3 for this gain was inserted into 
the gam matrix to obtain consistency with the matrix 
obtained for an angle of attack of 30° and to remove 
the necessity of scheduling gains while the airplane 
is stalled. The air data measurements necessary for 
the calculation of /? e st (eq (12)) probably cannot 
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be reliably obtained at extreme angles of attack 
Therefore, a lower limit of 50 m/sec was chosen for 
use in, the calculation For flight at airspeeds less 
than 50 m/sec, the following formula was used: 

a _ s 2m 

i ** " (JTTKfTI) Wpsc^ ( } 

The gain elimination process, discussed in this 
section, reduced the number of elements of the mea- 
surement vector feedback matrix from the original 18 
to the current 8. The cross coupling of the pilot con- 
trollers and the ideal rate response models was elimi- 
nated The feedback of the roll rate and the time in- 
tegral of the lateral acceleration were also removed 
The feedback of the turn coordination error to the 
directional pseudo control was found to be necessary 
at each of the angles of attack examined. This sig- 
nal was replaced by an estimate of the rate of change 
of the, sideslip angle at high angles of attack. Also 
at high angles of attack, the lateral acceleration was 
fed to the lateral pseudo control variable, and the 
yaw rate was fed to both pseudo control variables. 
The next step of the design process is to construct 
schedules which smoothly vary the gams as needed 
for changing flight condition. 

Gain Schedules 

In the preceding section, it was determined which 
gains generated by the solution of the LQR problem 
are necessary for acceptable stability and control of 
the airplane at angles of attack of 10°, 20°, 30°, 
and 60° In this section, this information is used 
in conjunction with the calculated gams given in 
figure 21 to construct gain schedules The schedules, 
along with the points examined in the preceding 
section, are given in figure 22. Those gains which 
were eliminated for all four angles of attack examined 
in the previous section are ignored m this section 
The remaining gains are scheduled as piecewise linear 
functions of angle of attack. The schedules generally 
follow the data of figure 21 with some deviations 
resulting from the gain elimination process of the 
preceding section. 

The calculated gams for the lateral stick deflec- 
tion (6 S ) and the roll rate model (p m ) to the lat- 
eral pseudo control variable (u] at ) are given in fig- 
ure 21(a). At low angles of attack (high airspeed), 
the gam on the roll rate model dominates. The pi- 
lot commands are filtered to limit the bandwidth of 
his stick motions which reach the controls At high 
angles of attack (low airspeed), the gain on the stick 
deflections dominates. Here the airplane is slow to 
respond to control, as is indicated by the small mag- 
nitudes of its eigenvalues. Therefore, the pilot stick 


motions are passed to the controls relatively unfil- 
tered. The schedules chosen for these gains are given 
in figure 22(a). At an angle of attack of 0°, the gain 
on the stick is 0.2, and the gam on the roll rate model 
is 0 8. At an angle of attack of 20°, the gain on the 
stick is 0 8, and the gain on the roll rate model is 
0 2. These gains are varied linearly at angles of at- 
tack between 0° and 20° and are constant at angles 
of attack beyond 20°. 

The calculated gams for the rudder pedal deflec- 
tion (<5 rp ) and the yaw rate model (r m ) are given in 
figure 21(b). These gains are nearly constant with 
changing angle of attack. Each has a value of ap- 
proximately 0.5. These gains and the yaw rate model 
dynamics were combined to form a lag-lead filter as 
follows: 

u dir = x _l i ^ r P (14) 

5 

This filter passes the low-frequency components of 
the rudder pedal commands (less than 5 rad/sec) di- 
rectly to the controls. The high-frequency compo- 
nents (greater than 10 rad/sec) are attenuated by 
6 dB 

The calculated gain from the turn coordination 
error (T ce ) to the directional pseudo control variable 
(udir) is given m figure 21(c). This gain varies 
between 1 2 and 1 8. A constant value of 1.4 was 
chosen As discussed in the preceding section, an 
estimated value of the rate of change of the sideslip 
angle (/? e st) was substituted for T ce for angles of 
attack of 20° or greater The transition between 
the two signals was specified as a linear function at 
angles of attack between 20° and 30°, as shown in 
figure 22(b). Simulation tests were conducted at an 
angle of attack of 15° by using (1) T ce only, (2) /? e st 
only, and (3) the scheduled mix of the two signals. 
No significant difference in the responses was noted 
among the three options. 

The calculated gains for the yaw rate (r) feed- 
backs are given in figure 21(c) Although these gains 
have values at all angles of attack, the gain elimi- 
nation process of the preceding section determined 
that they are effective only at the higher angles The 
gam schedules are given in figure 22(b). The sched- 
ules of the yaw rate gains are based primarily on 
the results of the preceding section and have little 
correlation with the calculated gains of figure 21(c). 
The scheduled gain from yaw rate (r) to directional 
pseudo control variable (u,j ir ) is linearly varied from 
0 to 0.30 over the angle-of-attack range of 20° to 30°. 
The yaw rate to lateral pseudo control variable (uj a t) 
is linearly varied from 0 to 0.60 over the angle-of- 
attack range of 10° to 30°. These gains are constant 
beyond 30°. 
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The calculated gains for the lateral acceleration 
( a y ) to the lateral pseudo control variable (u] at ) 
are given in figure 21(d). This gain increases from 
near 0 to approximately 0.18 over the angle-of-attack 
range of 10° to 30° and is fairly constant from 
30° to 60°. The schedule for this gam is given in 
figure 22(b). The schedule increases linearly from 
0 to 0.16 over the angle-of-attack range of 10° to 
25° and is constant beyond 25°. The gain for the 
lateral acceleration to the directional pseudo control 
variable (u<] ir ) is set to 0. 

Schematic diagrams of the lateral and directional 
control systems, including representations of the gam 
schedules and the control interconnect schedules, are 
presented in figure 23. A small-angle approxima- 
tion for angle of attack (a) has been applied to 
the formula for the turn coordination error (T C e) 
(eq. (2)) to produce the mechanization appearing in 
figure 23(b). Also, the value of velocity used in the 
T ce and 0 e st calculations has been limited to no less 
than 50 m/sec. 

Closed-Loop System Characteristics 

The eigenvalue characteristics of the airplane with 
the control system operating are plotted in figure 24. 
The absolute values of the eigenvalues and the damp- 
ing ratios of the second-order modes are presented as 
functions of angle of attack. These data may be com- 
pared with the data of figure 10, which are for the 
unaugmented airplane. 

The natural frequency of the Dutch roll mode 
is slightly less than 1.5 rad/sec, at angles of attack 
between 14° and 18° and less than 1.0 rad/sec above 
82°. The natural frequency has been increased near 
stall and in the poststall regime. For angles of 
attack near 90°, this increase is not present because 
low dynamic pressure prevents effective signals from 
the lateral accelerometer The damping ratio of 
the Dutch roll mode is maintained greater than 0.4, 
which is an improvement at angles of attack less than 
20° and greater than 50°. Between these regions, 
the damping ratio has been decreased. The damping 
factor (£w n ) has been increased by at least 70 percent 
for the data in figure 24. 

The coupled roll-spiral mode of the unaugmented 
airplane does not appear in figure 24. Instead, the 
spiral mode remains as a distinct first-order mode 
with a minimum time to half amplitude of 2.7 sec at 
an angle of attack of 30°. The roll mode couples with 
the differentiating filter used to generate /? e st The 
coupling occurs at an angle of attack of 24°, which is 
within the range where the turn coordination error is 
being replaced by Pest- The coupled mode exists up 
to 46°, and the damping ratio remains greater than 


0 95 For angles of attack greater than 46°, the first- 
order roll mode and the filter mode are again distinct. 
As the signal from the lateral accelerometer weakens 
with increasing angle of attack, the eigenvalue of the 
filter mode returns to its open-loop value. 

Airplane responses to pilot controller motions are 
presented in figure 25. The responses are for the 
linearized perturbation system with the control sys- 
tem operative. The stability axis roll rate responses 
to full-throw lateral stick deflection are given in fig- 
ure 25(a). The roll rates achieved are between 1.25 
and 1.60 rad/sec over the angle-of-attack range of 
10° to 70°. The responses for Jingles of attack of 20°, 
30°, and 40° have an initial delay of at most 0.30 sec. 
This is probably a consequence of the presence of 
the coupled roll-spiral mode m the unaugmented air- 
plane. At angles of attack of 10° and 70°, the de- 
lay is not present. In the former instance, the roll 
and spiral mode are distinct first-order modes with 
well-separated eigenvalues. The airplane is also being 
primarily driven through the roll rate model, which 
has a time constant of 1 sec In the latter instance, 
the aerodynamic forces and moments are ineffectual 
because of low dynamic pressure. The airplane pri- 
marily responds as a free body under the influence of 
forces and moments generated by the thrust-vector 
controls. 

The sideslip angle responses to full-throw rudder 
pedal deflections are given in figure 25(b). For the 
angle-of-attack range of 10° to 70°, the sideslip angle 
smoothly builds up with a small overshoot and no 
apparent oscillations. The sideslip angle achieved 
increases with increasing angle of attack. Lowering 
dynamic pressure causes smaller restoring moments, 
smaller aerodynamic moments, and smaller control 
moments because of smaller lateral accelerometer 
signals An initial delay of at most 0 3 sec is present 
in each response shown. For an angle of attack of 
10°, the sideslip angle decreases after peaking at 
2.5 sec because of the airplane rolling “away from 
the slip ” The roll rate can be arrested by applying a 
lateral stick deflection of approximately 50 percent 
in the “crossed-controls” sense. For the angle-of- 
attack range of 30° to 70°, the sideslip angle slowly 
builds up after the initial transient because of the 
airplane slowly rolling “into the slip.” Lateral stick 
motions of less than 7 percent in the “coordinated- 
controls” sense will arrest the roll rate in each of these 
instances. 

Concluding Remarks 

A preliminary design study of a lateral-directional 
control system using thrust-vectoring controls has 
been conducted. The control system design was 
based on a mathematical perturbation model of a 
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high-performance fighter airplane in trimmed flight 
over a wide range of angle of attack. 

Aerodynamic control surfaces and thrust- 
vectormg controls were scheduled so that they could 
be driven by separate lateral and directional pseudo 
controls This was accomplished by calculating com- 
binations of the controls which affected (1) the roll 
and spiral modes for the lateral pseudo control and 
(2) the Dutch roll mode for the directional pseudo 
control. During flight at low angles of attack and high 
airspeeds, aileron and differential tail movements are 
used for lateral control, and rudder movements are 
used for directional control. During flight at high 
angles of attack (near or beyond stall) and low air- 
speeds, rudder movements and yawing moments pro- 
duced by thrust vectoring are used for lateral control. 
Aileron and differential tail movements and rolling 
moments produced by thrust vectoring are used for 
directional control. 

State variable feedback gains were obtained from 
the solution of a linear quadratic regulator (LQR) 
problem which minimized the difference between air- 
plane angular rates and idealized rate models. The 
gams were used to calculate a set of sensor vari- 
able feedback gains. Of these gains, those which 
did not exhibit a beneficial effect upon simulated re- 
sponses to pilot control inputs were eliminated. The 
remaining gains were scheduled as functions of flight - 
condition 

The directional control channel uses a measure of 
the turn coordination error to improve the damping 
of the Dutch roll mode. During flight at high angles 
of attack, this feedback quantity is replaced by an 
estimate of the time rate of change of the sideslip 
angle, and yaw rate feedback is added The lateral 
control channel uses lateral acceleration and yaw 
rate during flight at high angles of attack Lag-lead 
prefilters are used on the control stick and rudder 
pedal signals. 

For the unaugmented airplane, the natural fre- 
quency of the Dutch roll mode is less than 1.0 rad/sec 
for angles of attack near 25° and greater than 77°. 
The damping ratio varies between 0 25 and 0 85. 
With the lateral-directional control system active, 
the augmented natural frequency is slightly less than 
1 5 rad/sec at angles of attack between 14° and 18° 
and is less than 1 0 rad/sec above 82°. The damping 
ratio stays greater than 0 40. The damping factor 


(£w n ) is increased by the control system at all angles 
of attack. 

For the unaugmented airplane, the roll and spiral 
modes merge to form a coupled roll-spiral mode for 
angles of attack greater than 18°. The damping ratio 
has a minimum value of 0 094 at an angle of attack of 
40° With the control system active, the augmented 
spiral mode remains a distinct first-order mode with 
a minimum time to half amplitude of 2 7 sec at 
an angle of attack of 30°. The roll mode interacts 
with the filter used to estimate the rate of change of 
sideslip angle (m the directional channel) to form a 
second-order mode at angles of attack between 24° 
and 46°. The minimum damping ratio of this coupled 
mode is 0.964. 

Thrust vectoring provides moments necessary for 
airplane control at extremely high angles of attack. 
The formation of pseudo control variables results in 
a system structure which is divided into lateral and 
directional components. Application of multivari- 
able controls design methodology produces feedback 
control loops which modify the airplane dynamic 
modes to make piloted control possible. Simulated 
responses to step pilot control inputs are stable and 
well behaved. For lateral stick deflections, peak sta- 
bility axis roll rates are consistently between 1.25 and 
1.60 rad/sec over an angle-of-attack range of 10° to 
70°. For rudder pedal deflections, the roll rates ac- 
companying the sideslip responses can be arrested by 
lateral stick motions of less than 7 percent. 


NASA Langley Research Center 
Hampton, VA 23665-5225 
July 11, 1985 
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TABLE I LINEAR PERTURBATION SYSTEM 
[a = 20°, V = 59 m/sec (115 knots)] 



r p i 


'—1.48 

0 

2.43 -5.92' 


r v 


- 144 

.254 

2 35 

614 

.656' 


o a 

d 

4> 


1.00 

0 

364 0 


<f> 


0 

0 

0 

0 

0 


0 r 

dt 

r 

" _ 

.010 

0 

- 449 - 017 

< 

r 

- + 

-.105 

-425 

039 

-.033 

-1 67 

< 

$D 

X 


1/9 J 


. .342 

156 

-.940 - 112. 


UJ 


. 001 

016 

-.005 

0 

.062. 


^v,roll 
A 


i, yaw 


a y>cg = -6 62 p + [ 030 501 - 104 0 954] { 


6a 

Sr 

6d 

,ro]l 


u 


u,yaw 


Variable 

Unit 

Range, 
rad (deg) 

6a 

rad 

±0 70(±40) 

6r 

rad 

±0.52(±30) 

6d 

rad 

±0.35(±20) 

6v,to\1 

rad 

±0 52(±30) 

6v, yaw 

rad 

±0.26(±15) 


Variable 

Unit 

P 

rad/sec 

<t> 

rad 

r 

rad/sec 

P 

rad 

a y, eg 

m/sec 2 
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TABLE II. MODIFIED PERTURBATION SYSTEM* 
[a = 20°; V = 59 m/sec (115 knots)] 


p 


'—1.48 

0 

2.43 

-5 92' 

<i> 


1 00 

0 

.364 

0 

r 

' 

.010 

0 

-.449 

-017 

UJ 


. 342 

156 

-.940 

- 112. 



■- 292 
0 

.650 

.-.024 



fly, eg — — 6 62/3 — 1.42t^ a t 052u ( jir 


f 1 


'-9.82 

33 O' 

6r 


-30 0 

-.024 

1 <5d 

. = 

6.00 

14.7 

I ^v,roll 


0 

19.0 

v ^v,yaw - 


.-13 9 

0 . 



*Units of measurement are given in table I 


TABLE III. LQR STATE VARIABLE FEEDBACK GAINS* 
[a = 20°, V = 59 m/sec (115 knots)] 







r p ) 







f Pm ) 



.308 

.128 

- 1.44 

. 785 ' 


\ 4> 1 

- + 

.110 

890 

-.213 

. 577 ' 

4 

1 <5s 1 

- + 

'. 0114 ' 

-.343 

-.181 

.777 

.139 

< 

1 r 

.341 

244 

.392 

.608 


I'm I 

.0083 













. <5 rp J 





and r m in rad/sec; y$ and V in m/sec; units for other variables are given in table I. 



TABLE IV. MEASUREMENT VARIABLE FEEDBACK GAINS* 


a, 

deg 

Pseudo 

control 

P, 

rad/sec 

t t 

J ce> 

rad/sec 

r, 

rad/sec 

^dy, 

m/sec 2 

?/3) 

m/sec 



r mi 

rad/sec 

<5rp 

10 

Uat 






0 38 





u dir 


-1.51 

KE9 







20 


B 

0.85 

§-0 78 

-0.118 


0.10 

0.90 

-0.22 

0 57 



Bill 

-1.16 

-0 33 

-0.021 


0.34 

0.25 

0 39 

0.61 

30 

*4at 

■Ui 

0 50 

-1 13 

-0 180 




US DEEM 



^dir 


-1.28 

-0 31 

0 001 





Wmm 

60 

«lat 


0 76 

-0 59 

-0 188 





0.45 


^dir 


-1.40 

#-0.02 

-0 025 




EH 

0.48 


*Gain elements remaining after the elimination process are given in boldface type, 
t Substitute 0 e st for T ce . Use equation (12) for a = 30° and equation (13) for a = 60°. 

■t-The lateral accelerometer is located 2.22 m forward of the center of gravity at the center of percussion. 
§Gain element was retained after elimination process, but substitute value of 0.3 was used instead of 0.78. 
# Value of 0 3 was used for consistency with gains for a = 30°. 
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Figure 1 Lift and drag coefficients. 6 ^ = 0°. 
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Figure 3 Horizontal-tail coefficients 
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Figure 6 Thrust-vector control moments. 
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Figure 7 Longitudinal trim. 
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Dutch roll mode — / 
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Figure 13 Aerodynamic control interconnects for lateral control. 
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Control Perturbation 



Figure 19 Modified perturbation system. 
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(b) Rudder 





Figure 21. Continued. 
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(a) Pilot controls and rate models. 
Figure 22. Gam schedules. 
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(b) Airplane motion feedback gains. 
Figure 22. Concluded. 


*>■ 



(a) Lateral channel. 
Figure 23. Control systems. 



All schedules are functions of angle 
of attack unless otherwise indicated. 
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Time, sec 

(a) Normalized stability axis roll rate response to lateral stick step input. 
Figure 25. Airplane responses to pilot controls. 


47 




1 Report No 

NASA TM-86425 

2 Government Accession No 

| 

3 Recipient’s Catalog No 

4 Title and Subtitle 

Preliminary Design Study of a Lateral-Directional Control System 

5 Report Date 
November 1985 



6 Performing Organization Code 

505-34-03-02 

7 Author(s) 


Frederick J Lallman 


8 Performing Organization Report No 

L-15990 

9 Performing Organization Name and Address 

NASA Langley Research Center 
Hampton, VA 23665-5225 

10 Work Unit No 

11 Contract or Grant No 

12 Sponsoring Agency Name and Address 
National Aeronautics and Space Administration 

13 Type of Report and Period Covered 
Technical Memorandum 

Washington, DC 20546-0001 


14 Sponsoring Agency Code 


15 Supplementary Notes 


16 Abstract 


A preliminary design of a lateral-directional control system for a fighter airplane capable of controlled 
operation at extreme angles of attack is developed. The subject airplane is representative of a modern 
twin-engine high-performance jet fighter. It is equipped with ailerons, rudder, and independent horizontal- 
tail surfaces. Idealized bidirectional thrust-vectoring engine nozzles are appended to the mathematic model 
of the airplane to provide additional control moments Optimal schedules for lateral and directional pseudo 
control variables are calculated. Use of pseudo controls results in coordinated operation of the aerodynamic 
and thrust-vectoring controls with minimum coupling between the lateral and directional airplane dynamics. 
Linear quadratic regulator designs are then used to specify a preliminary flight control system to improve 
the stability and response characteristics of the airplane Simulated responses to step pilot control inputs 
are stable and well behaved. For lateral stick deflections, peak stability axis roll rates are between 1.25 
and 1 60 rad/sec over an angle-of-attack range of 10° to 70°. For rudder pedal deflections, the roll rates 
accompanying the sideslip responses can be arrested by small lateral stick motions. 


17 Key Words (Suggested by Authors(s)) 18 Distribution Statement 

Airplane control Unclassified— Unlimited 

Integrated control 
Thrust-vectoring control 
Optimization 
Controllability 

Subject Category 08 

19 Security Classif (of this report) 20 Security Classif (of this page) 21 No of Pages 22 Price 

Unclassified Unclassified 49 A03 


For sale by the National Technical Information Service, Springfield, Virginia 22161 


NASA-Langley, 1985 























End of Document 



